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Abstract.  We  study  a  general  class  of  Euler  equations  driven  by  a  forcing  with  a  commu¬ 
tator  structure  of  the  form  [£,  u](p)  =  C(p\i)  —  C(p) u,  where  u  is  the  velocity  field  and  C  is 
the  “action”  which  belongs  to  a  rather  general  class  of  translation  invariant  operators.  Such 
systems  arise,  for  example,  as  the  hydrodynamic  description  of  velocity  alignment,  where 
action  involves  convolutions  with  bounded,  positive  influence  kernels,  =  (f>  *  f.  Our 

interest  lies  with  a  much  larger  class  of  £’s  which  are  neither  bounded  nor  positive. 

In  this  paper  we  develop  a  global  regularity  theory  in  the  one-dimensional  setting,  con¬ 
sidering  three  prototypical  sub-classes  of  actions.  We  prove  global  regularity  for  bounded  </>’ s 
which  otherwise  are  allowed  to  change  sign.  Here  we  derive  sharp  critical  thresholds  such 
that  sub-critical  initial  data  (po,uo)  give  rise  to  global  smooth  solutions.  Next,  we  study 
singular  actions  associated  with  C  =  —  (— dxx)a/2,  which  embed  the  fractional  Burgers’ 
equation  of  order  a.  We  prove  global  regularity  for  a  £  [1,  2).  Interestingly,  the  singularity 
of  the  fractional  kernel  \x\~<Jl+a\  avoids  an  initial  threshold  restriction.  Global  regularity  of 
the  critical  endpoint  a  =  1  follows  with  double-exponential  W1,00-bounds.  Finally,  for  the 
other  endpoint  a  =  2,  we  prove  the  global  regularity  of  the  Navier-Stokes  equations  with 
density-dependent  viscosity  associated  with  the  local  C  =  A. 
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1.  Fundamentals.  Euler  equations  with  a  commutator  structure 

We  are  concerned  with  a  new  class  of  Eulerian  dynamics  where  a  velocity  held,  u  :  x 
R+  i— y  Rn,  is  driven  by  the  system 


(1.1) 
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The  main  feature  here  is  the  commutator  structure  of  the  forcing 

(1.2)  T(pu)  =  [£,u](p)  :=  £(pu)  -£(p) u, 

expressed  in  terms  of  a  self-adjoint  operator  £  :  R  i— y  M  (the  action  on  pu  is  interpreted 
component-wise).  We  focus  on  the  Cauchy  problem  over  the  whole  space  =  Mn  or  over 
the  torus  =  Tn. 

A  typical  example  is  provided  by  radial  mollihers,  £(/)  =  (ft*  f,  associated  with  integrable 
(ft  G  L1,  which  yields  the  commutator  forcing 

(1.3)  T(p,u)(x)  =  (ft*  (pu)  -  (0*p)u  =  f  0(|x  —  y|)(u(y)  —  u(x))p(y)dy. 

J  i" 


The  corresponding  system  (1.1), (1.3)  arises  as  macroscopic  realization  of  the  Cucker-Smale 
agent-based  dynamics  [CS2007a,  CS2007b],  which  describes  the  collective  motion  of  N 
agents,  each  of  which  adjusts  its  velocity  to  a  weighted  average  of  velocities  of  its  neigh¬ 
bors  dictated  by  an  influence  function  (ft , 


Xj 


=  Vi, 

1  N 

=  ”xil)(vi  “  v*)’ 

i= i 


(xj,  Vj)  eTx  IT. 


For  large  crowds,  N  1,  one  is  led  to  the  hydrodynamic  description  (1.1), (1.3),  [HT2008, 
CCP2017].  For  recent  results  which  justify  the  passage  to  Cucker-Smale  kinetic  and  hydro- 
dynamic  descriptions  with  weakly  singular  kernels  <ft  (of  order  <  |)  we  refer  to  [Pes2015, 
PS2016].  The  global  regularity  of  such  one-  and  two-dimensional  systems  (1.1), (1.3)  was 
studied  in  [TT2014,  CCTT2016,  HT2016].  For  bounded,  positive  mollihers  it  was  shown  that 
there  exist  certain  critical  thresholds  in  the  phase  space  of  initial  configurations,  (po  >  0,  u0), 
such  that  sub-critical  initial  data  propagate  the  initial  smoothness  of  (p(-,  0),  u(-,  0))  = 
(po,  u0)  globally  in  time. 

Our  interest  lies  in  the  global  regularity  of  (1.1), (1.2)  for  a  much  larger  class  of  £’s  which 
are  neither  positive  nor  bounded.  We  have  three  typical  examples  in  mind. 


1.1.  Examples.  Consider  £  =  £^  of  the  form 

(!-4)  £*(/)(x)  :=  f  0(|x  -  y|)(/(y)  —  /(x))dy. 

J  R" 

Our  first  example  involves  bounded  kernels  with  a  finite  positive  mass,  denoted  (ft  G  Lj?  :  = 
{(ft  G  L°°  |  0  <  J  (ft(r)dr  <  cx)},  but  otherwise  are  allowed  to  change  sign.  The  resulting 
commutator  T<f>  =  [£^,u](p)  coincides  with  the  usual  convolution  action  in  (1.3), 

(1.5)  T(p,u)(x)  =  [£*,u](p)(x)  =  f  <K|x-y|)(u(y)-u(x))p(y)dy,  (ft  e  L%. 

J  R" 

The  action  £^  in  (1.4)  and  its  commutator  forcing  (1.5)  are  well  defined  for  non-integrable 
0’s  as  well.  As  a  second  example  we  consider,  0Q(x)  :=  |x|_(-n+a\  associated  with  the  action 
of  the  fractional  Laplacian1  £«(/)  =  —A a(f),  a  <  2, 

A a(/)(x)  =  p.v.  I  fjfjdy,  Aa  =  (-A)"/2,  0  <  a  <  2. 

jRn  lx  y| 

xWe  shall  abuse  notations  by  abbreviating  Cc/)a  :=  Ca  since  the  distinction  is  clear  from  the  context  of 
the  sub-index  involved. 
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The  corresponding  forcing  is  then  given  by  the  singular  integral 

(1.6)  T(p,  u)(x)  =  -Aa(pu)  +  Aa(p)u  =  p.v.  f  — ^^p(y)^y- 

jRn  ix  yi 

The  operator  T  in  (1.6)  is  well-dehned  as  a  distribution  over  the  whole  space  =  Mn.  When 
dealing  with  the  torus  =  Tn,  the  forcing  T  can  be  expressed  in  terms  of  the  periodized 
kernel  <j>a{ z)  =  £keZri  |z+2Jk|n+a. 

Finally,  as  a  third  example  we  consider  the  full  Laplacian  C  —  A  corresponding  to  the 
limiting  case  a  =  2  with  forcing  7~(p,  u)  =  pAu  +  2 (Vp  •  V)u.  This  leads  to  the  density- 
dependent  system  of  pressureless  Navier-Stokes  equations 


(1.7) 


Pt  +  V  •  (pu)  =  0, 

(pu)t  +  V(pu  (8)  u)  =  V  (p2Du), 


Du  =  {diUj}. 


We  close  by  noting  that  these  equations  are  typically  come  “equipped”  with  certain  standard 
global  bounds.  Thus,  in  addition  to  the  obvious  conservation  of  mass, 


M0:= 


J  p0(x)dx  = 


p(x,  £)dx, 


we  have,  since  C  is  assumed  self-adjoint,  J  (p£(pu)  —  £(p)pu)dx  =  0,  conservation  of  mo¬ 
mentum,  JRn  pu(-,  t)  dx  =  fRn  p 0u0  dx.  Also  for  C $  we  have  the  p-weighted  energy-enstrophy 
bound 

p(x)p(y)0(|x-y|)|u(x)-u(y)|2dxdydt  =  /  p0|u0|2dx. 

xl"  it" 


(1.8)  f  p|u|2dx+  [  f 
d]R™xIT|  Jo  J] Rn 


1.2.  The  one  dimensional  case.  Statement  of  main  results.  The  main  focus  of  this 
paper  is  one-dimensional  case  where  (1.1), (1.2)  reads, 


(1.9) 


Pt  +  (p«) 

(pu)t  +  (pu2) 


0, 

pC(pu)  -p£(p)u, 


(. x ,  t)  G  x  M+. 


We  shall  make  a  detailed  study  on  the  propagation  of  regularity  of  (1.9)  for  sub-critical 
initial  data,  dictated  by  the  properties  of  C. 

We  begin  by  recalling  that  (1.9)  with  C  =  amounts  to  the  one- dimensional  Cucker- 
Smale  “flocking  hydrodynamics”  [CS2007a,  CS2007b,  MT2014,  CCP2017] 


(1.10) 


Pt  +  (pu)x  =  0, 

(pu)t  +  (pu2)x  =  /  (j)(\x -y\)(u(y)  -  u(x))p(x)p(y)dy, 


(. x ,  t)  e  D  x  M+. 


Global  regularity  for  bounded  positive  0’s  persists  if  and  only  if  the  initial  data  are  sub-critical 
in  the  sense  that  [CCTT2016] 

(1.11)  u'Q(x)  +  (j)*  po(x)  ^  0  for  all  i6R. 

In  section  3.1  we  extend  this  regularity  result  for  general  bounded  0’s  whether  positive  or 
not.  The  results  below  are  stated  over  the  torus,  f!  =  T1,  for  the  purely  technical  reason  of 
securing  a  uniform  lower  bound  of  the  density  away  from  vacuum,  which  in  turn  provides 
uniform  parabolicity  of  the  w-equation.  However,  the  local  well-posedness  follows  from  our 
analysis  over  =  M  line  as  well. 
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Theorem  1.1.  Consider  the  hydrodynamics  flocking  model  (1.10)  with  a  bounded  mollifier, 
(ft  G  Lfft  having  a  positive  total  mass  I  {(f)  =  f  <ft{r)dr  >  0,  and  subject  to  sub-critical  initial 
data  (po,Uo)  G  (L^T1),  W/1,0O(T1))7  such  that 

u'0(x )  +  (ft*  p0{x )  >0,  x  G  T1. 


Then  (1.10)  admits  global  smooth  solution. 

Next,  we  extend  this  result  to  the  case  of  the  positive  singular  modifiers  (fta(r)  =  |r|”d+Q) 
Pt  +  {pu)x  =  0, 


(1.12)  ^  2  f  u{y)  —  u{x)  .  .  .  . 

'  (pu)t  +  {pu  )x  =  p.v.  I  — — „.i  1+J P(x)p{y)dy,  a  <2, 


(. x ,  t)  G  12  x  M+. 


\x  -  y | 


Here  we  follow  a  general  iteration  scheme  for  proving  (higher)  regularity  outlined  in  section 
2,  in  which  one  seeks  bounds  on  the  density,  p,  and  then  bounds  the  “action”  C{p).  The 
uniform  bounds  on  the  density  for  all  three  cases  are  worked  out  in  section  3.  We  then  turn 
to  secure  bounds  on  the  action  or  —  what  amounts  to  the  same  thing,  uniform  bound  on  ux, 
which  in  turn  yields  global  well-posedness.  In  section  4  we  discuss  the  global  regularity  for 
bounded  modifiers  C  =  in  (1.10),  and  in  section  5  for  the  the  Navier-Stokes  equations, 

t  =  a, 


(1.13) 


Pt  +  {pu)3 
{pu)t  +  {pu2)a 


0, 

{p~ ^x)  XI 


(x,  t)  G  fl  x  M+. 


This  is  the  one-dimensional  special  case  of  the  general  class  of  Navier-Stokes  equations  stud¬ 
ied  in  [BGDV2007], 

Theorem  1.2.  Consider  the  Navier-Stokes  equations  (1.13)  subject  to  initial  data  {uq,  pf)  G 
H2{ T1)  x  H3{ T1).  Then  (1.13)  admits  a  global  solution  in  the  same  class. 

Finally,  in  section  6  we  prove  the  global  smooth  solutions  for  the  commutator  forcing 
associated  with  the  singular  action  Ca{p )  =  —A a{p)  corresponding  to  singular  kernel  (fta{r )  = 
|r|~(1+Q),  1  ^  a  <  2. 

Theorem  1.3.  Consider  the  system  of  equations  (1.12)  with  1  ^  a  <  2  subject  to  initial 
data  (uo,  po)  G  H3( T1)  x  H2+a{ T1).  Then  (1.12)  admits  a  global  solution  in  the  same  class. 

It  is  remarkable  that  the  singularity  of  (fta  =  |x|_^1+Q')  removes  the  requirement  for  a  finite 
critical  threshold  which  is  otherwise  called  for  integrable  (ft  G  Lff.  Specifically,  in  section 
3.2  we  prove  that  for  any  singular  kernel  such  that  limr^o[mod27r]  r  ■  min|z|<r.  (ft{\z\)  f  oo,  the 
density  of  the  corresponding  system  (1.10)  remains  uniformly  bounded  which  in  turn  drives 
the  global  regularity.  The  analysis  of  equation  with  the  singular  action  Ca  becomes  critical 
when  a  reaches  value  1.  The  necessary  kF1,00-bounds  on  the  solution  pair  {u,  p)  in  this  case 
admit  double-exponential  growth  in  time,  consult  (6.22). 
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2.  Propagation  of  global  regularity.  A  general  iteration  scheme. 


2.1.  L°°-bound  of  the  velocity.  We  assume  that  £  satisfies  the  following  monotonicity 
condition.  Let  x+  =  arg  max  g  (x)  and  =  argmin^(x).  Then  for  /  ^  0 


(2.1) 


f  C(fg)(x+)  <  £{f)(x+)g(x+),  g(x+)  =  maxg(x) 

I  £{fg){xJ)  ^  £{f){x-)g{x-),  g{xJ)  =  minr/(x) 

V  X 


which  holds  for  C  =  with  positive  0’s.  Application  of  (2.1)  with  (f,g)  =  ( p,u )  implies 
that 

„  argrnax  w(-,  f) 

T(p,u)(x+)  <  0  ^  T{p,u){x-),  x±  = 


arg  min  u(-,t) 

and  yields  that  u  in  (1.9)  (and  likewise  —  the  velocity  components  Ui  in  (1.1))  satisfy 
maximum / minimum  principle 

(2.2)  minuj(x,  0)  ^  Mj(x,  t)  ^  maxuj(x,  0) 

X  X 

Likewise,  ||u(-,£)||x,°o  remains  finite  for  0  G  L #  with  arbitrary  sign. 

2.2.  Critical  threshold  and  a  first  order  conservation  law.  We  outline  our  general 
strategy  for  tracing  the  global  regularity  of  (1.9).  The  key  observation  is  that  the  commutator 
form  of  (1.9)  entails  the  transport  of  ux  +  £(p)  away  from  vacuum.  To  this  end,  differentiate 
(1.1)  to  find  that  u'  :=  ux  satisfies 

(2.3)  u't  +  uu'x  +  ( v! )2  =  £(pu)x  —  u£(p)x  —  u'£(p). 

For  the  latter  we  use  the  density  equation,  £(pu)x  =  £((pu)x)  =  —£(p)t  to  conclude 

{v!  +  £{p))t  +  u(v!  +  £(p))x  +  u(u  +  £{p))  =  0. 

This  calls  for  introduction  of  the  new  variable,  e  :=  u'  +  £(p),  which  is  found  to  satisfy 

(2.4)  et  +  (ue)x  =  0,  e  =  u'  +  £(p). 

Together  with  the  density  equation,  this  yields  that  e/p  is  governed  by  the  transport  equation 


(2.5) 


P 


-  +«  -  =0. 


P 


Hence  e/p  remains  constant  along  the  characteristics  x(t)  =  u(x(t),t), 

(2  e(x(t),t)  =  e0(x) 

p(x(t),t)  Po(x)  ' 

It  follows  that  if  eo/po  is  allowed  to  have  singularities,  then  these  initial  singularities  will 
propagate  along  characteristics  and  a  solution  of  (1.9)  will  consist  of  strips  of  regularity 
trapped  between  the  curves  carrying  these  singularities.  To  avoid  this  scenario,  calls  for  the 
following  bound  to  hold. 

Assumption  2.1.  [Critical  threshold]  There  exist  finite  constants  p_  ^  0  <  r/+  such  that 

e0(x) 


(2.7) 


V- 


< 


Po{x) 


^  p+  for  all 
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Remark  2.2.  We  note  in  passing  that  integration  of  (2.7)  yields  rj-M0  f  (V0  +  C(po))dx. 
Hence,  since  £<j,(po)  has  zero  mean  and  u(-,  t )  is  either  periodic  or  assumed  to  have  vanishing 
far-held  boundary  values,  it  follows  that  (2.7)  requires  r/_  ^  0. 

We  will  investigate  the  propagation  of  regularity  of  solutions  subject  to  sub-critical  initial 
data  (2.7). 

2.3.  The  iteration  scheme  —  a  priori  control  estimates  via  e.  The  study  of  global 
well-posedness  for  all  three  cases  of  commutator  forcing  we  have  in  mind  —  bounded,  sin¬ 
gular  and  local  (NS)  mollihers,  share  a  common  scheme  of  establishing  control  over  the  key 
quantities,  even  though  the  handling  of  the  three  cases  is  quite  different  when  it  comes  to 
analytic  details.  In  this  section  we  highlight  those  main  common  features  in  three  steps. 

•  Step  #1  ( Pointwise  bounds  on  the  density).  Our  aim  is  to  show  that  for  a  certain  range 
of  threshold  bounds  p_  ^  0  <  ?/+,  the  density  remains  bounded  from  above  and  away  from 
the  vacuum 

(2.8)  0  <  p-  ^  p(-,f)  ^  p+  <  oo. 

In  view  of  transportation  of  the  ratio  e/p,  (2.6),  we  also  have 

(2.9)  ^  ^  *7-  <  °- 

We  conclude  that  the  quantity  of  interest,  e  =  ux  +  £(p),  will  remain  uniformly  bounded, 

e_  :=  p_p+  ^  e(-,t)  ^  e+  :=  p+p+. 

•  Step  #2  ( Pointwise  bound  on  the  action  £(p)  and  slope  ux).  Equipped  with  the  uniform 
bound  on  e  we  turn  to  establish  a  bound  on  the  action  £(p),  which  is  equivalent  to  controlling 
the  slope  ux.  In  the  case  of  bounded  mollihers  we  seek  a  pointwise  bound  on  the  action  £(p) 

(2.10)  £_^£(p)^£+,  pe  L°°. 

This  will  imply  the  desired  CT-bound  of  the  velocity 

V-P+  -  £_  C  ux(-,t)  £  ?/+p+  +  £+. 

For  singular  fractional  mollihers  £a,  we  focus  on  the  critical  case  a  —  1,  where  we  use  a 
nonlocal  maximum  principle  to  establish  control  over  p'  which  in  turn  enables  us  to  control 
ux  indirectly,  thus  avoiding  an  additional  obstacle  coming  from  the  Hilbert  transform.  For 
the  NS  case,  we  hrst  control  the  slope  ux  via  energy  bounds,  then  conclude  with  control  of 

£(p)  Pxx- 

ft  is  clear  from  the  fact  that  the  higher-order  quantity  e  satishes  lower-order  estimates 
that  a  proper  statement  of  well-posendess  result  for  singular  mollihers  requires  p  to  be  in  a 
regularity  class  A"s+Q  provided  u  is  in  the  class  A"s+1,  while  e  is  in  the  class  Xs. 

•  Step  ^3  ( Higher  regularity  control).  The  necessary  bounds  sought  in  (2. 8), (2. 10)  may 
require  a  restricted  set  of  initial  configurations  depending  on  finite  critical  threshold  assumed 
in  (2.7).  Whether  these  thresholds  )]±  are  restricted  or  not,  the  corresponding  bounds  will 
be  derived  solely  on  the  basis  of  the  mass  equation  for  p,  and  the  fact  that  e  =  ux  +  £(p) 
satishes  the  transport  equation  (2.5).  This  argument  can  be  iterated  to  higher  derivatives 
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as  follows.  Note  that  if  a  quantity  Q  is  transported,  Qt  +  uQx  =  0,  then  the  same  transport 
equation  governs  Qx/ p 


(2.11) 


=  0. 


Let  us  apply  this  argument  to  Q  =  e/p:  then  if  \{e/p)x\/p  is  bounded  at  t  —  0  it  will  remain 
bounded  at  later  time.  Unraveling  the  formulas,  we  obtain  the  pointwise  bound 


(2.12) 


e'(x,t)\  <  C(e±,p±)\p'(x,t)\. 


This  control  bound  will  become  a  key  tool  in  proving  Theorem  f.3. 

Now  that  (e/p)x/ p  is  transported,  we  can  apply  the  argument  above  repeatedly  to  obtain 
a  hierarchy  of  pointwise  bounds 

(2.13)  |e^(x,t)|  ^  C\p(k\x,  £)|,  k  —  0,1,... 

It  is  therefore  clear  that  such  bounds  would  allow  to  apply  the  same  control  principle  as 
stated  above  in  extending  our  results  into  higher  order  Sobolev  spaces.  However,  we  will 
leave  to  pursue  this  direction  to  a  future  work. 


3.  Bounded  density  in  one-dimensional  equations  in  commutator  form 

In  this  section  we  implement  the  above  strategy  for  global  regularity  in  the  presence  of 
commutator  forcing,  70,  depending  on  the  properties  of  the  modifier  0.  We  begin  with  a 
general  discussion  on  the  boundedness  of  the  density  sought  in  step  #1.  Here,  the  bound 
(2.8)  is  driven  by  the  diffusive  character  of  the  mass  equation,  which  is  revealed  once  we 
rewrite  the  mass  equation  of  (1.9)  in  the  form 

(3.1)  Pt  +  upx  =  -ep  +  p£(p). 

In  view  of  the  uniform  bound  (2.9),  we  see  that  ep  ~  p2  behaves  as  a  quadratic  term.  This 
implies 

(3.2)  -  //+p2  +  p£(p)  ^  pt  +  upx  ^  -p_p2  +  p£(p). 

We  turn  to  check  step  #1  in  the  three  cases  of  interest. 


3.1.  Bounded  density  with  bounded  modifiers  £^,  0  E  L#.  Consider  the  case  of  C  = 
Cfj,  —  f  0( \x  —  y\)(p(y)  —  p(x))dy  with  0  E  L which  is  assumed  to  have  a  positive  mass 
f  (ft(r)dr  >  0.  We  emphasize  that  0  need  not  be  positive.  We  verify  the  boundedness  of  p 
using  the  straightforward  bound 

(3.3)  -  J(0)p  -  |0|ooMo  ^  £^(p)  ^  -!{4>)p  +  HooMa,  /(0)  :=  J  0(r)dr  >  0. 

Inserted  into  (3.2)  we  find 

-(?/+  +  J(0))p2  -  |0|oo M0p  ^  Pt  +  Upx  ^  -(p_  +  J(0))p2  +  |0|oo M0p. 

The  inequality  on  the  left  shows  that  along  characteristics,  the  density  is  bounded  away  from 
vacuum  by  a  lower-bound  p(t )  >  e^^°°Motpo.  The  inequality  on  the  right  shows  that  the 
density  remains  bounded  from  above  for  any  p_  >  — J(0),  that  is,  provided  (2.7)  holds  for 
such  p_’s, 


(3.4) 


u'Q[x)  +  0  *  po(x)  -  /(0)po(x)  ^  P-Po(^), 


V-  >  -!{<!>)■ 
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3.2.  Bounded  density  with  singular  mollifiers  Ca,  a  <  2.  To  bound  the  density  from 
above,  we  consider  the  case  of  positive  mollifiers  which  are  singular  in  the  sense  that 


(3.5)  lim r 777,0 (r )  t  cc,  m^r)  min0(|z|). 

In  this  case  we  use  the  bound 


£(p)(x+)<  <f>(\x  -  y\)(p(y)  -  p+)dy 

J  \x—y\^r 

^  7770 (r)  /  (p(y)  -  p+)dy  ^  ui^Mq  -  2rm(f)(r)p+. 

J \x— 

By  assumption,  for  any  7/_  ^  0  we  can  choose  a  small  enough  r  =  r+  such  that  2r+7?70(r+)  = 
1  —  p_  and  the  bound  on  the  right  of  (3.2)  then  implies  that  the  maximal  value  of  the  density 
p+(t)  =  p(x+(t),t )  satisfies 

P+  <  ~V-P2+  -  2r+7770(r+)p+  +  c0p+  ^  -p+  +  c0p+,  c0  =  m0(r+)Mo 


Thus,  p(-,t)  remains  bounded  from  above.  We  conclude  that  for  singular  kernels  satisfy¬ 
ing  (3.5)  ,  the  density  remains  upper-bounded  independent  of  the  lower  threshold  r/_ .  In 
particular,  this  applies  to  (j>a(r)  =  r~^1+a\a  <  2. 

We  turn  to  the  lower  bound  on  the  density  away  from  vacuum.  For  positive  0’s,  whether 
singular  or  not,  we  have; 

C4>(p)(x-)=  [ 0(\x  ~  y\)(p(y)  -  p(x-))p(y)dy  ^  0. 
dy 

Therefore,  the  inequality  on  the  left  of  (3.2)  implies  that  minima  values  of  the  density, 
p~(t)  =  p(x-(t),t)  at  any  interior  point  X-(t)  =  a,rgmm^R{p(y,t)}  with  |x_|  <  R,  satisfy 
p_  ^  — P+p2.  and  hence  p(-,£)  >  0.  In  the  particular  case  of  the  torus  fl  =  T1,  we  conclude 
with  a  uniform  lower  bound  away  from  vacuum 

(3.6)  p(-,t)  >  p-(t)  =  - — ..  ,  (po)_  =  minpo(x)  >  0. 

tp+(p0)_  +  1 


3.3.  Bounded  density  with  NS  equations  £2  =  dxx.  We  use  the  regularization  coming 
from  the  parabolic  part  of  the  mass  equation  which  becomes  evident  when  (1.13)  is  written 
in  the  form 


(3.7)  pt  +  upx  +  ep  =  ppxx. 

It  implies  the  lower-bound  which  could  be  read  from  the  LHS  of  (3.2),  p_  ^  —p+p2_  which 
recovers  the  same  lower-bound  (3.6).  Trying  to  pursue  the  same  argument  for  an  upper- 
bound  of  the  density  fails  when  using  the  RHS  of  (3.2).  Instead,  we  note  that  the  quantity 
/  :=  u  +  px  is  the  primitive  of  e  and  hence  satisfies  the  transport  equation  ft  +  ufx  =  0. 
This  follows  by  direct  computation  of  (1.9) 

ut  +  uux  —  -(p't  +  up'x ). 


It  follows  that 


(3.8)  |p'(-,f)|  ^  2|t70|oo  +  IPoloo 


2This  is  a  special  case  of  the  monotonicity  condition  (2.1)  with  (f,g)  =  (1 ,  p)  implies  C(p)( X-)  ^ 
C(l[x-))p-  =  0. 
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Here  and  throughout  |  •  |p,  1  ^  p  ^  oo,  denotes  the  Lp-norm.  Given  the  uniform  bound  on 
p'  and  since  we  already  proved  that  p  >  0,  (3.8)  ties  the  upper  bound  for  p  as  well.  We  note 
in  passing  that  even  though  we  can  now  express  the  density  equation  as  a  pure  diffusion 

(3.9)  pt  =  ppxx  +  F, 

with  bounded  forcing  F  =  —  upx  —  pe  G  L°°,  we  can  only  reach  the  end-point  Schauder 
estimate  pxx  G  BMO  (see  [Schll996]),  which  is  not  enough  to  secure  a  uniform  bound  of 
£(p)  =  pxx  necessary  to  get  control  over  the  slope  ux.  We  will  provide  additional  details 
how  to  reach  that  bound,  which  is  needed  for  the  global  existence  of  NS  equations  in  section 
5  below. 


4.  Global  existence:  bounded  mollifiers,  £  = 

With  regard  to  Theorem  1.1,  it  is  straightforward  to  verify  step  in  the  case  of  bounded 
mollifiers  —  in  view  of  (3.3),  the  upper-bound  of  p  implies  that  £(p)  is  uniformly  bounded, 

| £^(p)  +  /(0)p+ 1  ^  MocMo,  and  hence  ux  =  e  —  £(p)  is  uniformly  bounded.  We  conclude 
the  global  regularity  for  sub-critical  initial  data  satisfying  (3.4),  namely,  for  a  fixed  e  >  0 
there  holds 

u'Q(x)  +  (p*  p0(x )  ^  ep0(x),  e  >  0. 

In  the  particular  case  of  T1,  this  requires  the  positivity  of  u'0  +  f>*  p0  stated  in  theorem  1.1. 
This  recovers  the  same  critical  threshold  of  positive  mollifiers  (1.11). 

5.  Global  existence:  Navier-Stokes  equations,  £  =  A 

In  this  section  we  will  prove  Theorem  1.2.  Recall  that  the  boundedness  of  p  conveys  to 
boundedness  of  e  =  ux  +  pxx,  via 

(5.1)  \e(x,t)\  ^rj\p(x,t)\, 

and  that  p  satisfies  further  a  priori  Cl -regularity  (3.8).  Note  that  these  low-regularity  a 
priori  bounds  hold  classically  under  the  assumptions  u  G  H2,  p  G  H3,  which  are  the  spaces 
for  which  Theorem  1.2  is  stated,  and  these  are  the  lowest  integer  regularity  //"-spaces  that 
justify  the  above  computations.  We  now  proceed  by  establishing  a  priori  estimates  in  these 
spaces. 

First,  let  us  quantify  control  over  the  high-order  regularity  of  e. 

Lemma  5.1.  For  each  n  —  0,1, ...  we  have  the  following  a  priori  estimate 

(5.2)  dt\e\ 2Hn  ^  C( \e\2Hn  +  M|n+i)(|w*|oo  +  |e|oo)- 

Proof.  For  n  =  0  the  Lemma  follows  easily  by  testing  the  e  equation  (2.4).  For  n  —  1, ...,  let 
us  differentiate  (2.4)  n  times  and  test  with  e^n\  We  obtain  (dropping  the  integral  signs) 

n+1 

dt\e{n)\l  <  we(n+1)e(n)  +  ^«(fc)e(n+1_fc)e(n). 

k= 1 

For  the  first  term  we  integrate  by  parts  to  obtain  trivially  |ue^"+1')e^n')|  <  liGloclel^n.  For 
each  of  the  remaining  terms  on  the  right,  k  =  1, . . .  ,n  +  1,  we  apply  Gagliardo-Nirenberg 
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inequalities,  | dlf\*n  ^  |/| 


n  \  +  \  n 


I  w 


.,  1  ^  i  ^  n,  obtaining 


(fe)„(n+l-fe)  (n)| 


it  'e 


^  le^Uw^U  \e 


(n+1— fc)  | 


2n 


^  \e{n)\2\u 


1- 

rr  loo 


le(n) 


2n+l  —  fc 


fc- 1 
fc-1  fc-1 

n  L,  I  n 

I  UX  |  }{  n 

fc-1 


,1- 
^  I  oo 


n+1  — fc 

n+1  — fc  n+1— fc 

n  lei 


\Ux\Hn  \Ux\oo 


I  I1” 

e  oo 


Hn 

n+1  — fc 


and  by  Young’s  inequality  l+UgG+i  ^  (\e\2Hn  +  \u\2Hn+1) 

the  proof. 


\u3 


,+  e 


which  completed 

□ 


We  now  proceed  establishing  bounds  on  ux  and  uxx  in  a  sequence  of  increasing  norms, 
which  eventually  will  close  the  estimates  together  with  Lemma  5.1.  Recall 

(5.3)  ut  +  uu'  =  pu"  +  2  p'u' . 


Testing  with  u  and  using  (3.8)  we  obtain 


dt\u\\  = 


p\u'\  +  /  p'uv!  ^  --  /  p\u\2  +  2p-_1\u\l 


This  proves  the  natural  energy  bound  u  G  D  L2H^.  Next,  we  test  with  —  u"  to  obtain 

(dropping  the  integrals) 


dt\u'\2  =  \uu'u"\  —  p\u"\z  +  2|  p'u'u 


< 


I  // 1 2  i  I  1 2  I  /|2  i  I  / 1 2 1  / 1 2 

-p\u  |  +  \u\  \u  |  +  | p  |  | u  I  . 


Using  uniform  bound  on  u  and  (3.8), 


dt\u\l  ^  — -p_ P-  =  min  p{'i  t )  >  0, 

which  implies  u  G  L^H],  H  L2H2.  In  particular,  this  implies  |mx(-+)|oo  £  L1,  and  hence 
the  estimates  on  H 1  -norm  of  e  from  Lemma  5.1  closes  with  an  integrable  multiplier  on  the 
right  hand  side  of  (5.2).  It  remains  to  establish  a  further  similar  bound  on  |m"|2  to  close  the 
estimate  on  the  grand  quantity  \u"\\  +  l^~KII  +  l  p\2hs-  So,  we  differentiate  (5.3)  twice 
and  test  with  u" : 

a  // 1 2  ,  J  /  //  //  /  //N  /  III  ,  r,  III  I  II  I  A  II  II  II  ,  n  I  III  H 

Ot\u  |2  +  -u  u  u  =  — ( pu  )  u  +  2puu+Apuu+2puu 

I  /// 1 2  I  I  III  II  I  O  II I  I  II  i  A  II  II  II 

=  —  p\u  |  +  pu  u  +2 p  uu  +  4p  u  u  . 


So, 


dt\u"\l  < 


\ux 


\u 


// 1 2 


-p-  1 11 


III  1 2 


+  1  p' 


oo 


u 

n\2 


II 1 2 
2 


+  I  |  oo  |  p//7 1 2 1  ^ 1 2  +  4|e 


un 


\u 


n  1 2 


P-W"\i  +  kU|e+  +  K+) 


Given  the  established  integrability  of  Im^Ioo  and  Lemma  5.1  we  have  proved  boundedness  in 
H 2  for  u,  and  H3  for  p. 
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6.  Global  existence:  singular  mollifiers  C  =  CQ,  1  ^  a  <  2 

In  this  section  we  prove  global  regularity  result  for  the  equation  with  fractional  C  =  —  Aa  in 
space  of  data  H 3  x  H2+a.  The  case  a  =  1  is  critical  similar  to  the  classical  fractional  Burgers 
equation,  [KNS2008,  CV2010,  CV2012]  but  with  additional  non-linearity  in  the  dissipation 
term.  We  will  leave  the  subcritical  case  a  >  1  as  an  easy  consequence  of  the  proof  presented 
here  for  the  case  a  —  1.  Note  that  with  the  initial  datum  (uq,  po)  in  H 3  we  can  avoid  making 
assumptions  on  eo  as  eo  G  H 2  C  C1+1  for  any  7  <  1/2  by  the  Sobolev  embedding.  With  this 
we  recall  a  priori  uniform  bounds  from  the  previous  section, 

(6.1)  sup  |e|oo  <00,  0  <  p_  ^  p(x,t )  ^  p. |_,  w_  ^  u(x,t )  ^  u+. 

0  <t<T 

on  any  finite  time  interval  of  existence.  The  lower  bound  on  the  density  is  the  main  reason 
why  we  resort  to  the  periodic  domain.  In  the  open  space  such  bound  is  only  known  to  hold 
on  any  finite  interval,  lacking  a  uniform  parabolicity  to  the  system. 

Moreover,  for  a  solution  in  H 3  the  transport  equation  (2.11)  for  Q  =  e/p  can  be  solved 
classically  along  characteristics  of  u  which  results  in  the  bound  (2.12)  which  we  quote  for 
convenience 


(6.2) 


ex(x,t)\  ^  C\px(x,t)\,  for  all  (x,t)  G  [0,T)  x  T. 


The  proof  will  consist  of  four  steps.  First,  we  establish  local  existence  in  H3  by  obtaining 
rough  a  priori  bounds  without  exploiting  dissipation  term.  This  allows  to  perform  classical 
desingularization  of  the  kernel  as  an  approximate  scheme  to  obtain  local  solutions.  Second, 
we  establish  uniform  control  over  first  order  quantities  (pxloo,  l^zloo  over  the  interval  of 
regularity.  The  strategy  here  resembles  the  treatment  of  the  critical  SQG  by  Constantin 
and  Vicol  [CV2012],  but  with  additional  technicalities  related  to  the  non-linear  nature  of 
the  dissipation  term.  We  then  invoke  the  results  of  Schwab  and  Silverstre  [SS2012]  to  obtain 
instantaneous  C7-regularization  and  use  it  to  have  an  easier  control  on  the  oscillations  in  the 
midrange  of  scales  of  the  non-linearity.  Third,  we  establish  uniform  control  over  H2  norm 
of  solutions  by  proving  an  analogue  of  the  Beale-Kato-Majda  estimates.  With  the  H 2  and 
IT1’00  bounds  we  finally  conclude  by  a  proving  a  uniform  control  of  the  penultimate  H3- norm 
of  the  solution  on  the  entire  interval  of  existence. 

It  will  be  useful  to  introduce  the  following  notation.  For  three  functions  /,  g,  h  of  x ,  z  we 
denote 


Hf,g,h)  :=l 


f(x,z)g(x,z)h(x,z)  dzdx 

Izl2 


Moreover,  for  a  cutoff  function  p  and  parameter  r  >  0  we  denote 


®<r(f,9,h)  =  ^ 
®>r{f,9,h)  =  ^ 


f(x,z)g{x,z)h(x,z) 


p(z/r)dzdx 


f(x,z)g(x,z)h(x,z) 


(1  —  p(z/r))dzdx. 


In  the  sequel  we  will  also  use  the  following  notation  8zf(x )  =  f(x  +  z)  —  f(x),  and  the 
expansion 


(1  -8)f"(x  +  9z)d9. 


(6.3) 


SJ(x)  =  f'(x)z  +  z 2 


0 
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6.1.  Local  well-posedness  in  H 3:  a  priori  estimates  without  the  use  of  dissipation. 

The  purpose  of  this  section  is  to  obtain  a  priori  estimates  in  H3  which  do  not  rely  on  the 
dissipation  term.  Namely,  we  will  obtain  the  classical  Riccati  equation  for  the  quantity 

1  =  \u\h3  +  \p\h3  ~  M H3  +  Ie| H2  +  | pW- 

Yt  ^  CY 2, 

which  is  independent  of  desingularization  of  the  kernel  Kx  =  - - — s+r.  This  allows  to 

(|2|2+52)t~ 

conclude  local  existence  via  the  classical  approximation  methods. 

Let  us  write  the  equation  for  u'": 

(6.4)  <  +  uvl”  +  Au'u'"  +  3 u"u"  =  T{p'",  u )  +  3 T{p",  u')  +  37V,  «")  +  T{p,  um). 
Testing  with  v!"  we  obtain  (we  suppress  integral  signs  and  note  that  J  u"u"u"'  =  0) 

(6.5)  dt\u'"\l  =  -7 u\u"')2  +  2 T{p'",  u)um  +  6 T{p",  u')um  +  67V,  u")K"  +  2^(P> 

We  will  now  perform  several  estimates  with  the  purpose  of  extracting  term  u'"  \  \  on  the  right 
hand  side,  times  a  lower  order  term  in  u  and  possibly  a  top  order  term  in  p  which  we  will 
address  subsequently.  First,  we  have  trivially 

(6.6)  \u\u"f\  ^  In'IoV"!!- 
Let  us  estimate  the  dissipative  term  first: 

j  T(p,u'")u"'dx  =  r)(U"'M-«"'(U)r^4Y 

Switching  x  and  y  and  adding  cross-terms  p(x)u"'(x)  we  obtain 

jnP,u'")u"'dx=-1-  J [  p(x)(u'"(y)  - 

+  \jJ (pW  -  P(xW\y)  -  u'"(x))^L. 

The  first  term  is  clearly  negative.  We  note  in  passing  that  it  is  bounded  below  by 

Jf  P{x){u'"{y)  -  u'"(x)f  ^y_^2  ^  P-W"\2hi/2. 

While  it  is  undoubtedly  a  crucial  piece  of  information,  it  does  depend  on  the  fact  that  the 
kernel  is  singular.  As  we  indicated  earlier,  however,  we  seek  estimates  that  are  independent 
of  singularity.  So,  at  this  point  we  will  simply  dismiss  the  dissipation  term.  As  to  the  second 
term,  we  rewrite  it  as 

JJ u"(x)(p(y)  -  p(x)){u"' (y)  -u"\x))-^-^  =  $(u"',5zp,5zu"'). 

We  estimate  the  large-scale  part  of  the  integral  using  integrability  of  \z\~2  at  infinity  as 
follows 

T>1(u///,4p,4yA//)  ^  \p\ooW\l- 


(6.7) 
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As  to  the  small  scale,  we  use  the  expansion  (6.3)  on  p.  We  have 
$<i (u'",8zp,8zu'")  =  [[  y(z)u"\x)p'(x)5zu"\x)d"  dx 


+  J  (1  —  9)  J  J  ip(z)u'" (x) p" (x  +  9z)8zu"'[x)dz  dx  d6. 
Writing  the  first  integral  in  the  principal  value  sense  results  in  the  cancellation 


cp(z)u'"(x)p'(x)u'"(x)- —  =  0, 


while 


tp(z)u'"(x)p'(x)u"(x  +  z)d~  dl  =  [  u,"(x)p'(x)Hip(u'")(x)dx, 


where  H $  is  the  truncated  Hilbert  transform  given  by  convolution  with  the  kernel  .  It 
a  bounded  operator  on  L2.  We  thus  have  the  estimate 


u"’  (x)pr  (x)Hv  (■ u ){x)dx 


Putting  the  estimates  together,  we  arrive  at  the  bound 


(6.8) 


T{pyy"dx^c\u'X(\p\oo  +  \p'\oo) 


We  now  proceed  with  the  remaining  three  terms  in  (6.5)  in  a  similar  fashion.  We  have 

T(p"',  u)u'"  =  $(/'(•  +  z), 8zu )  =  $>i(pw(-  +  z),  um,  Szu)  +  $<i (p'"(-  +  z),  5zu) 
^  W"\2\p"'\2\u\oo  +  J  Hip(pm)(x)u'"(x)u'(x)dx 

+  [  (1  —  9)  f  f  p"\x  +  z)u'"{x)u"  (x  +  9z)tp(z)dz  dx  d9 


<  Klal/S'laduloo  +  Moo  +  KI*). 


T(p",  u')v!"  =  $(/(•  +  z),  u"\  5zu')  =  $>!(/(•  +  z),um,  8zu')  +  /  Hy)(x)u'"{x)u'\x)d 


+  J  (1  —  9)  j  j  p"(x  +  z)v!" (x)u" (x  +  9z)ip(z)  dz  dx  d9 
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And  the  last  term  requires  more  preparation, 
T{p'y)u"'=  ffp'{y)u"\x)(u"{y)-u"{x)) 


dy  dx 


x  -  y\‘ 


(/ o'{y)u'"{x )  -  p\x)u"\y))(u'\y)  -  u'\x)) 


dy  dx 

\x  -  y I2 


{p\y)-p\x))u"\x){u"(y)-u'\x)) 

p\x){u"\x)-u"\y)){u"{y)-u"(x)) 


dy  dx 

\x  -  y\ 2 

dy  dx 

\x  -  y I2 


=  ^<3>(<5zp,,u,,,,<5zu")  -  ,8zv!" ,8zu") 


dz 


=  yySzp  ,u",Szu”)  +  -  /  /  <p(z)p'(x)u "  (; x)5zu  (x)—  dx 


+  -  /  (1  —  9)  /  /  ip(z)p  (x  +  9z)u  (x)5zu  (x)dz  dx  d6 


d  z 

p'(x)({s,y(x)fy  jjj  dx 


<  |p'|0o|«'"|2|«"|2  +  |/|=o|«'"|2|«"|2  +  |/'|2|«"'|2|«"|»  +  7  /  /  P» 


The  latter  integral  is  bounded  by  |p"|oo|u"|^.1/2  ^  |p"|oo|u/"|2|u//|0o-  Putting  the  obtained 
estimates  together  we  obtain 

$|«'"ll  <  ckiKIu'U  +  IpIoo  +  Ip'Ioo  +  \p"U 

(6.9)  +  K'|2(|p"|2  +  |  p"  1 2  +  |p"|oo)(Moo  +  K|oo  +  Kloo) 

+  |u///|2|u"|2(|p,|oo  +  Ip'Ioo)- 


Finally,  by  Sobolev  embedding,  |u"|2  +  |u|oo  +  (w'loo  +  |u"|oo  ^  C{\u"'\2  +  |u|2),  and  |p|oo  + 
Ip'Ioo  +  | p" 1 00  +  | p" 1 2  <  C(\p'"\2  +  |p|2)  which  results  in  the  bound 

(6.10)  dt\u"'\l  <  (K |2  +  \u\2)2{\p"'\2  +  |p|2)  +  (K|2  +  M2)3. 

To  control  the  energy  \u\2  we  avoid  using  the  natural  balance  relation  (1.8).  Instead  we  test 
(1.1)  directly  with  u.  Performing  much  the  same  estimates  as  above  we  obtain,  for  example, 

dt\u\\  ^  \u\oo\u\2\p\2  +  Ip'hKhMoo- 

Putting  this  together  with  (6.10)  we  obtain  the  Riccati  equation  for  the  H3- norm: 

(6.11)  dt\u\H3  <  \u\H3\p\H3  +  \u\2h3. 

In  order  to  close  the  estimates  we  now  have  to  find  a  similar  bound  on  the  H3- norm  of  p. 
This  cannot  be  done  directly  by  manipulating  with  the  density  transport  equation.  Instead 
we  will  make  use  of  the  transport  of  the  first  order  quantity  e,  in  terms  of  which  we  will 
provide  the  final  estimates.  Let  us  note  the  inequality 

\p\h3  ^  \u\H3  +  \e\H2  +  | p 1 2 - 


dt\u\H3  ^  \u\h3(\^\h2  +  \pU)  +  \u\%3- 


Thus, 

(6.12) 
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From  Lemma  5.1  we  have  the  bound  on  |e|#2: 

(6.13)  <9t|e|#2  C(\e\H2  +  |u|#3)2. 

And  the  similar  bound  holds  for  |p|2.  We  have  obtained  the  classical  Riccati  equation  for 
the  quantity  Y  =  \u\H3  +  \e\H2  +  |p|2: 

Y,  Y  CY2. 

Note  that  Y  ~  \u\h 3  +  \p\h3,  hence  we  have  proved  necessary  a  priori  bound  for  the  local 
well-posedness  in  H3. 

6.2.  Control  over  (n^loo  and  \px\oo  on  intervals  of  regularity.  Suppose  that  we  have  a 
classical  solution  (u,  p)  G  C([0,T);f/3)  as  proved  to  exist  in  the  previous  section.  We  now 
seek  to  establish  a  uniform  bound  on  lu^oo  and  \px\oo  on  the  entire  interval  [0,T).  First, 
let  us  recall  that  we  have  already  established  a  priori  uniform  bounds  of  e,p  and  u  in  terms 
of  the  finite  initial  quantities  e±,p±  and  u±,  consult  (6.1).  Next,  as  we  noted  the  density  p 
satisfies  a  parabolic  form  of  the  density  equation: 

(6.14)  pt  +  upx  +  ep  =  pC(p) 

Similarly,  one  can  write  the  equation  for  the  momentum  m  =  pu : 

(6.15)  mt  +  umx  +  em  —  pC(m) 

Note  that  in  both  cases  the  drift  u  and  the  forcing  ep  or  em  are  bounded  a  priori.  Moreover, 
the  diffusion  operator  has  kernel 

K{x,h,t)  =  p(x)jf^ 

which  satisfies  all  the  assumptions  of  Schwab  and  Silverstre  [SS2012],  A  direct  application 
of  [SS2012]  tells  us  that  there  exists  an  7  >  0  such  that 

|p|c,T'(Tx[T/2,T))  ^  C(\p\l°°(0,T)  +  \pe\L°°(0,T)) 

(6.16)  [m\c-y(Tx[T/2,T))  ^  +  |^e|ioo(o,T)) 

|m|ct(Tx[T/2,T))  ^  C{\u\l°°(P,T)i  |p|l°°(0,T))) 

where  the  latter  follows  from  the  first  two  since  p  is  bounded  below.  Of  course,  since  u,  p  are 
in  H3  on  [0,  T)  this  implies  C'7- bound  on  the  entire  interval  of  regularity,  however  we  need 
the  bound  to  be  independent  of  H3,  which  may  blow  up,  in  the  second  half  of  it.  It  is  also 
interesting  to  note  that  the  original  equation  for  u  has  a  kernel  K(x,  h,  t )  =  p(x  +  h)  not 
even  with  respect  to  h,  so  no  known  results  on  regularization  are  directly  applicable  to  the 
w-equation. 

Remark  6.1.  In  regard  to  higher  order  regularization  via  Schauder,  we  make  the  following 
observation.  For  Q  =  e/p  we  recall  that  Qx  was  shown  to  be  under  control  (note  that  this 
still  doesn’t  imply  that  either  ex  or  px  are  under  control).  Hence,  trivially,  \Q\c-y  remains 
bounded  at  all  times.  Denote 

jf  rw  \  Q(x  +  h)  -  Q(x) 
fihQ(x)  =  - — - , 

and  note 

A  n(v\  =  She(x)  4-  e(x)5hp(x) 
h  p{x  +  h)  p{x  +  h)p{x) 
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Since  p  is  C 7  and  bounded  away  from  zero  this  implies  that  e  G  C1  with  |e(t)|c"?  ^  C/t1 . 
With  this  in  mind,  we  now  have  the  momentum  equation  in  the  form 

mt  +  b(x)mx  +  a(x)Am  =  F, 

where  the  drift  b,  the  coefficient  function  a  and  the  source  F  are  all  in  C1 .  As  of  this  writing 
there  has  been  no  known  Schauder-type  bounds  proved  for  an  equation  in  such  generality 
despite  many  recent  developments  in  that  cover  partial  cases,  see  [CK2015,  DZ2015,  SS2012, 
TJ2015].  The  question  presents  an  independent  interest  and  we  will  address  it  in  subsequent 
work. 


Let  us  now  establish  control  over  p' .  We  write 

dtp'  +  up"  +  up  +  e'p  +  ep  =  —p'Ap  —  pAp' . 

Using  again  u'  =  e  +  A p  we  rewrite 

dtp'  +  up"  +  e'  p  +  2  ep'  =  —2p'Ap  —  pAp' . 

Let  us  evaluate  it  at  the  maximum  of  p'  and  multiply  by  p'  again  (we  use  the  classical 
Rademacher  theorem  here  to  justify  the  time  derivative): 

(6.17)  dt\p'\2  +  e!  pp  +  2e|p'|-  =  -2|p'|2Ap  -  pp’Ap. 

In  view  of  (6.1)  and  (6.2)  we  can  bound 

| e' pp  +  2e\p'\2\  ^  Cjp'|2. 

Next,  using  the  nonlinear  bounds  from  [CV2012]  we  have 

(6.18) 
where 

Dp\x)  = 


PP'Ap'  ^  \p-Dp\x)  +  cP-\p'\l  ^  CiDp'(x)  +  c2|p'|L 
4  p+ 


| p'{x)  -  p'(x  +  z)\‘ 


dz. 


Using  smooth  decompositions  of  the  underlying  R  in  all  of  the  below  we  have 


A  p(x)  =  Ftp  = 


f  p'(x  +  z)  -  p'(x) 

I \z\<r  % 

f  p(x  +  z)  —  p(x) 

’2tt<\z\  \  - 


dz 


dz. 


p(x  +  z)  —  p(x) 


dz 


'  r<\z\<2iT 


The  latter  is  clearly  bounded  by  a  constant  C3  depending  only  on  p_,  which  in  (6.17)  results 
simply  in  the  bound  c3|p,|^0.  The  first  is  bounded,  via  Holder,  by 


\p'\loVrD1/2p'(x)  ^  -c1Dp,(x)  +  c4r\p,\ 


4 

00  ’ 


Note  that  this  term  gets  absorbed  by  the  dissipation  (6.18)  entirely  if 

c2 


r  = 


4c4|p'| 


The  integral  in  the  middle  is  bounded  by,  using  U7-regularity, 

Ip'lLlPlc-r/r1-7  =  C5|p'|^7  <  c6  +  flp'lL, 
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where  the  cubic  term  again  is  absorbed  by  the  dissipation.  Putting  the  estimates  together 
we  obtain 

(6.19)  dt\p'\2  ^  c6  +  c3|p'|-  -  c7Dp'(x), 

which  establishes  the  claimed  control  of  p' .  We  intentionally  keep  the  dissipation  term  as  it 
still  will  be  used  on  the  next  step  to  absorb  other  terms. 

Now  we  can  do  the  same  for  the  momentum  derivative  mx.  Clearly  it  is  sufficient  to  finish 
the  proof  for  ux  as  well.  Note  that  the  equation  for  momentum  is  similar,  so  we  will  skip 
details  that  are  similar.  We  have 

dtm'  +  um"  +  u'm!  +  e'rri  +  em'  =  —p'  Am  —  pAm' . 

Evaluating  at  maximum,  multiplying  by  rn\  and  using  bounds  on  e,  e!  we  have 

(6.20)  dt\m'\2  ^  c8(|m'|^  +  Iplool^'U)  +  \'m'\2\Ap\  +  \p'\\m'\\Am\  -  c9Dm'(x )  -  Ci0|m'|^. 

As  to  \p'\ Im'IIAm)  we  proceed  as  before,  loosing  p'  in  view  of  already  established  control  over 
it.  We  obtain  the  bound  simply  by  taking  r  =  1: 

Cu\m'\  +  \m,\D1^2rm' (x)  ^  Cu\m'\  +  Ci2|m,|2  +  Dm!(x ), 

with  the  latter  being  absorbed  again  in  the  dissipation.  As  to  the  term  | /rz'  | J  | Ap |  we  still 
proceed  as  before,  however  in  the  mid-range  integral  r  <  \z\  <  2ir  we  use  the  full  force  of 
the  obtained  bound  on  p'.  This  results  in  logarithmic  optimization  bound 

|m'|2|Ap|  ^  Ci3|m'|2(l  +  lnr  +  y/rD1^2 p' (x)). 

Ignoring  the  trivial  quadratic  term  \m!  |2  we  have 

cislm'l2  lnr  +  ci3\m'\2y/rD1/2p\x)  ^  cislm'l2  lnr  +  cu\m'\4r  +  Dp'(x ). 

Notice  that  the  latter  will  be  absorbed  by  the  dissipation  term  in  (6.19)  when  we  add  the 
two  equations  together.  Choosing 


2cu\m'\ 

we  obtain  for  the  In  r  and  r-terms  above  the  bound 

Ci5|m  |  In  \m  \  +  ~^\m  I  5 

with  the  latter  being  absorbed  into  the  cubic  term  in  (6.20).  Altogether  we  have 

(6.21)  dt\m'\2  ^  ci6|m'|2(l  +  ln+  \m'\)  +  ^ Dp\x ). 

We  now  have  to  add  the  two  equations  (6.21)  and  (6.19)  together  to  absorb  the  residual 
Dp'- term  and  obtain  the  final  bound 

(6.22)  dt{\m'\2  +  |p'|2)  ^  c17(\mf  +  \p'\2)(l  +  ln+(|mf  +  |pf)). 

This  implies  double-exponential,  but  finite,  bound  on  the  given  interval.  This  also  finishes 
the  proof. 
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6.3.  Control  over  H 2  via  kF1,0°.  In  this  section  we  will  establish  an  estimate  on  the  H2- 
norm  of  the  solution 

X  ~  F  |2  +  Ip  I2  F  I2  +  le  I2 

in  terms  of  \ux\oo  is  a  manner  similar  to  the  Beale-Kato-Majda  criterion.  Namely,  we  will 
prove 

(6.23)  X'  ^  C{  1  +  |'R,|oo)X(l  +  log+  X). 

Given  the  result  of  the  previous  section  this  establishes  uniform  bound  in  H2  on  the  interval 
of  existence  [0,T)  of  an  if3-solution.  The  equation  for  u"  reads 

u"  +  mm"  +  3 u'u"  =  T(p",  u )  +  2 T(p\  u ')  +  T(p,  u"). 

Testing  with  u"  the  local  terms  ,  after  integration  by  parts  ,  become  bounded  by  X  | v!  \  ^  triv¬ 
ially.  We  now  look  into  key  estimates  for  the  right  hand  side.  We  will  start  with  what  proved 
to  be  the  most  involved  term  in  the  previous  section.  We  skip  the  standard  symmetrization 
and  addition  of  cross-product  terms  in  the  calculations  below  and  typically  display  the  final 
representations.  We  have 

(6.24)  /  /  T{p\u')u”dxdy  =  &(8zp' ,  8zu' ,u")  +  ^(p',8zu',8zu''). 


For  the  second  term  we  have  8zu'8zu"  =  l((8zu')2)x.  So,  switching  the  derivative  onto  p '  we 
obtain 

<I>{p,  8zu,  8zu")  =  — ^*h(p",  8zu',  8zu). 

Now,  we  bound  small  and  large  scale  parts  as  follows 

\$>r(p",8zu',8zu')\  ^  -\p"\2\u'\l, 

r 

and 

\^<r(p",8zu,8zu)\  ^  Vr\p"\2\u\lV3/4,4, 

where  in  the  latter  we  used  the  Holder  and  Gagliardo-Sobolevskii  definition  of  IF3/4,4  space. 
Optimizing  over  r  we  obtain 

\§{p' ,8zu' ,8zu")\  ^  |p"|2|u/|4/3|u/|^3/4.4, 
and  by  Gagliardo-Nirenberg, 

W\w^a  <  \nXlM2, 


and  interpolation  we  obtain 


,8zv! ,8zu")\  ^  |p,,|2|M/|2/3|u,|^3|u/,|i/i3/2  ^  -\p"\22\u'UW\oo  +  £ 


u 


// 1 2 


l/f1/2' 


With  £  <  P-1‘1  the  last  term  is  absorbed  by  the  dissipation.  Finally,  by  Gagliardo-Nirenberg 
we  have 


(6.25)  K|4  ^  K'l2/2ML/2- 

Recalling  that  |w|^i,oo  is  under  control,  we  finally  obtain 

\$(p',8zu',8zu")\  ^  C\p"\l/2\u"\l/2\u'\0O  +e\u”\2Hl/2  ^  CX  +  e\u"\Hi/2. 
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For  the  other  term  <&(8zp' ,  8zu' ,u")  the  splitting  is  necessary  but  optimization  is  not.  We 
have,  in  view  of  (6.25), 

$>1  (8zp',8zu',u")  ^  \p\i\u'\i\u"\2  <  \p"\12/2\u"\l/2  ^  X. 

As  to  $<i,  we  write  8zu'{x )  =  8~u'{x )  —  zu"(x)  +  zu"{x),  and  note  that  1 8zu'(x)  —  zu"{x)\  ^ 
\z\^2D1/2u"{x).  So,  we  have 


H \pp'  (x)\u"  (x)\2  dx 


+  |P'|oo|^"|2K|2 


^\H^u%  +  -\Prjuy2  +  e\u%l/2. 

Note  that  since  Hv  is  a  bounded  Fourier  multiplier,  by  the  log-Sobolev  inequality,  we  have 


So, 


\Hvfl\oo  <  |p/|oo(1  +  log+  \p"h)  <  |p'|oo(l  +  log+X). 

I <  1  ( 8 z p/ ,  8zi i' ,  u")\  ^  67(1  +  | p  1 00  +  \p'\lo  +  \u'UX(l  +  log+  X)  +  e\v!  \2hi/2- 


Taking  into  account  the  control  over  W1,0°  norms,  we  have  proved  the  bound 

1 7V, u')u"\  <  CX{  1  +  log_|_ X)  +  e\u"\Hi/2. 

Next,  let  us  bound  the  dissipation  term 

T(p,  u")u"  dydx  =  — <h(p,  8zu" ,  5/u")  +  <f>(<52p,  <52w",  u"). 


Obviously, 


As  to  &(8zp,8zu,,,u ")  we  have 


§(p,8zu",8zu")  >  p_|u"|^1/2 


$>r(8zP,  8ZU  ,U  )  ^  —\u  |2, 
r 


and 


Optimizing  we  obtain 


/  |«"(*)|  / 

/K  «/|z|<2r 

2 

dzdx  ^  V^’|p,|oo|m"|2|w"|j?i/2. 


s:  Ip^VlfV'lhh  <  +  (vi 

«e|«"iy/1+CA'(l  +  log+A'), 

which  closes  the  estimates  with  the  help  of  dissipation. 

It  remains  to  estimate  the  last  term.  By  switching  x  and  y  we  obtain 

T(p",u)u"  dydx  =  [[  p"(x)(u{x)  -  u(y))u"(y)-  dydX 


l?y"  1 2 

OO  |  W  |2 


\x-y  |2 

X(x)(u(x)  -«(»))(«"(»)  -«"(!))■  ** 


I®  -  2/|s 


F  _  2/ 1 

=  <&{p" ,  52u,  52m")  +  [  p"u"A(u)dx. 


(6.26) 
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Clearly,  by  the  log-Sobolev  inequality, 

|  f  p"u"A(u)dx\  ^  Ip'WhlMoc  ^  h'loo^l  +  log+X). 


For  the  F-term  we  have 


while 


\$>r(p",5zu,5zu")\  ^  -\p"\2\u"\2, 


\^<r(p",SzU,5zu”)\  ^  |m'|oo  J  | P"{%)\ 
Optimizing,  we  get 


\6zu"(x)\ 


\z\<2r 


dzdx  ^  \u'\00\/r\p"\2\u"\Hi/2. 


ny)\  <  |/|2|ii"|^/3|m'|^3|m"||?;2  <  e|u%,/s  +  -e\p"\TwtfW\oo  <  cx. 


We  have  proved  that 


dt\u”\i  ^  ~£\u"\2h1/2  +  CX(  1  +  log+X). 


As  to  quantity  e,  we  apply  Lemma  5.1  to  obtain 

dt\e’\l  sC  CX. 

Putting  the  estimates  together,  (6.23)  follows. 


6.4.  Control  over  H 3  via  H 2  and  W1,oc.  For  a  given  classical  solution  (w,  p )  E  C([0,  T);  i/3) 
we  have  established  uniform  bounds  on  \ux,  px |oo  and  |w,  p|#2  on  the  entire  interval  [0,  T ).  We 
now  seek  to  establish  final  control  over  the  iL3-norms.  Note  that  we  already  have  estimate 
(6.13)  which  with  the  new  information  readily  implies 


d+ 


< 


+  u" 


2 

2’ 


Now  we  get  to  bounds  on  \u'"\\.  Not  surprisingly  all  of  the  estimates  mimic  the  already 
obtained  sharper  estimates  for  H2  with  the  use  of  dissipation.  In  what  follows  we  will 
indicate  necessary  changes  and  refer  to  appropriate  places  in  Section  6.3  for  details.  Also, 
we  will  drop  from  the  estimates  all  quantities  that  are  already  known  to  be  bounded,  such 
as  \u,p\H2,  etc.  Thus,  following  (6.4)  we  can  see  that  all  the  terms  on  the  left  hand  side  obey 
the  bound  by  |w,|0o|m/,,|2  ^5  |u///||.  We  are  left  with  the  four  terms  on  the  right  hand  side: 


T{p'",u)u'",  T(p",  u')u'",  T(p' ,u")u'",  T(p,  u'")u"'. 


First,  the  dissipation  term  obeys  the  same  bound  (6.8)  where  we  now  keep  the  dissipation  : 


(6.27)  /  T(p,  u"')u"'dx  ^  -p_ 


u 


1 2 

I  H1/2 


T  C  |  v!"  1 2  ( |  P  |  oo  T  \p'\ 


< 


-p-\u 


1 2 

I  if1/2 


+  C\u" 


Next,  the  term  T{p'" ,u)u'"  will  be  estimated  in  the  same  way  as  (6.26)  with  replacements 
p"  — >•  p'",  u"  — >  v!" .  We  have  the  bound 


|T(p">K"|  ^  \p"'\2\u"'\2\Au\ 


+  10 


u 


1 2 

I H1/2 


10, 

+  — I PJ 
p- 


lll  1 3/2 


U 


III  1 1/2 


Since  (Auloo  ^  \u\h2  <  C  and  | p/// 1 2  ^  | e" 1 2  +  W"\2  we  have 


\T{p'"i  u)u"'\  < 


I2  + 

12  w 


U 


I2  4- 
12  w 


P_ 

10 


\u 


2 

m/ 2‘ 
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Next,  the  term  T(p" will  also  be  estimates  as  in  (6.26)  with  a  simple  replacement 
u  — *  u ',  i.e.  raising  the  derivative  of  u  by  one  on  every  step.  We  obtain  directly, 

l/T-/  II  l\  III  I  ^  I  III  I  lll\  I  A  /  I  I  I  III  I  2  |  //  1 3/2  I  ///  1  1/2  I  //I 

<  p  2  «  2  Ait  |oo  +  e\u  bri/2  +  -  p  2  u  I2  loo- 

e 

Dropping  |p"|2  and  using  that  lAit'l,*,,  |'U//|00  ^  |w|h3,  we  obtain 


iz-j-/  //  /\  ///1  ^  P-  1  /// 1 2  1 

|T(p,m)m  I  ^  —  |w  |^i/2  + 


it 


H3 


P-  1  /// 1 2 


r  —  \  rnyz  ,  /~i  ,  \  m\2 

T7T  lr/1/2  +  C  +  |m  |2. 


10 


Finally,  the  term  T(p7 ,u")u"'  can  be  estimates  as  term  (6.24)  by  raising  the  derivative  of  u 
by  one  and  with  the  use  of  boundedness  of  |p"|2,  |//|oo.  We  obtain 


1  <-rt  J  I  ^  „L.«/|2  I  i  //|3/2i  //i  I  //i  .  I  , 

I  /  (p , u  ju  1  ^  1^1/2  “t-  |p  I2  l^i  ol^  1 00  4“  |p 


w  1 1/2 1  m  1 3/2 


it 


u,,(x)Hp'(x)u'"(x)dx 


,1/1  I  , 

_  |p  |oo |n 


/// 1 2 

2- 


We  have  trivially,  |it//|4|-u//|0o  ^  |it|^3,  and 


u"(x)Hp'  (x)u'"  (x)dx 


^\u"\2\u’"\2\Hp'\00<\u'"\2\Hp"\2<\u'"\2. 


This  completes  the  estimate  for  the  H3- norm  Y:  Y'  ^  CD  on  the  time  interval  of  existence. 
This  completes  the  proof. 
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